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Abstract 

The conditions for the existence of force-free non-relativistic translationally invariant one- 
dimensional (ID) Vlasov-Maxwell (VM) equilibria are investigated using general properties of 
the ID VM equilibrium problem. As has been shown before, the ID VM equilibrium equations 
are equivalent to the motion of a pseudo-particle in a conservative pseudo-potential, with the 
pseudo-potential being proportional to one of the diagonal components of the plasma pressure 
tensor. The basic equations are here derived in a different way to previous work. Based on this 
theoretical framework, a necessary condition on the pseudo-potential (plasma pressure) to allow 
for force-free ID VM equilibria is formulated. It is shown that linear force-free ID VM solutions, 
which so far are the only force-free ID VM solutions known, correspond to the case where the 
pseudo-potential is an attractive central potential. A general class of distribution functions leading 
to central pseudo-potentials is discussed. 

PACS numbers: 52.20.-j, 52.25.Xz, 52.55.-s, 52.65.Ff 
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I. INTRODUCTION 



Plasma equilibria are suitable starting points for investigations of, for example, plasma 
instabilities and plasma waves. For collisionless plasmas, the relevant equilibria are self- 
consistent solutions of the stationary Vlasov-Maxwell (VM) equations (see e.g. Refs. 1,^). 

In the present paper we focus exclusively on non-relativistic one-dimensional quasi-neutral 
VM equilibria with translational symmetry, with the distribution functions depending only 
on the Hamiltonian and the two canonical momenta corresponding to the invariant di- 
rections (in this paper chosen to be the x- and y-directions). A large amount of work 
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especially on one-dimensional current sheets and plasma boundary layers, which are of fun- 
damental importance for the structure and stability of plasmas as many plasma activity 
processes, e.g. magnetic reconnection,— happen there preferentially. 

In the present paper it is not our aim to add to the plethora of ID VM solutions existing 
already, but to use some generic properties of the ID VM equilibrium problem to investigate 
the conditions for the existence of force-free ID VM equilibria (see e.g. Ref. 1361 ) . An 
obvious property of ID VM equilibria is that the structure needs to be in force balance. 
In the quasi-neutral case this means that the sum of the magnetic pressure and one of the 
diagonal components of the pressure tensor has to be constant. The relevant component 
of the pressure tensor is the one for the single spatial coordinate upon which the equilibria 
depend. In this paper we will choose this coordinate to be z, so the component of the 
pressure tensor in the force balance equation will be P zz . In a number of papers it has also 
been noticed that the translationally invariant ID VM problem is equivalent to the motion 
of a pseudo-particle in a conservative pseudo-potential and/or that the force balance for the 
ID VM structure is equivalent to pseudo-energy conservation (see e.g. Refs. 
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29c the pseudo-particle analogy has also recently been used for MHD 



Directly connected to the pseudo-particle analogy and the related pseudo-energy conser- 
vation law (equivalent to force balance of the ID VM equilbrium) is the special role played 
by P zz - As was first noticed by Grad 4 for the case of vanishing electric potential and only 
one non- vanishing magnetic field and vector potential component, but otherwise arbitrary 
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distribution functions, the derivative of P zz with respect to the non-vanishing component 
of the vector potential equals (modulo constants) the current density. This was generalized 
by Bertotti^ who included a non-vanishing electric potential and showed that the partial 
derivative of P^ with respect to the electric potential is proportional to the charge density 



fsee also Refs. 



371 for the same conclusion for 2D VM equilibria) . Lerche^ then generalized 



this to equilibria with two magnetic field, vector potential and current density components, 
but for a restricted class of distribution functions. Channel!^ investigated the case of van- 
ishing electric potential and a special class of distribution functions for which he showed that 
the two components of the current density are, again modulo constant factors, given by the 
partial derivatives of the particle density with respect to the components of the vector po- 
tential. If investigated carefully one can see that the full expressions including the constant 
factors are again the partial derivatives of P zz with respect to the components of the vector 
potential. Using the force balance condition, Mynick et al.— then showed that, indepen- 
dently of the distribution function, the partial derivatives of P zz with respect to the electric 
potential and the two components of the vector potential are always proportional to the 
charge density and the components of the current density, respectively. They also showed 
that this property is maintained under the assumption of quasi-neutrality. The authors 
then use this property to construct ID VM equilibria with certain defined properties. The 
special role of P zz in the context of ID VM equilibria has also been emphasized by Attico 
and Pegoraro^, again for the case of vanishing electric potential. Similarly to Channel!^, 
they used this property to construct a number of special distribution functions for ID VM 
equilibria. More recently, Mottez 3 -^ gave a detailed discussion of the role of the full pressure 
tensor (not just of one component) for distribution functions of the same type as discussed 
by Channel! 21 , but including the case of non- vanishing electric potential. 

The properties of P zz make it a very useful quantity to start any investigation of ID 
VM equilibria, since all other quantities such as particle density, charge density and current 
density can be derived by differentiation. P zz is also equivalent to the pseudo-potential of 
the analogous pseudo-particle problem and can thus be used to study the properties of ID 
VM equilibria qualitatively without the need to solve the equilibrium differential equations. 
In particular, the pseudo-particle analogy can be used to formulate conditions on P zz that it 
has to satisfy to allow the existence of ID force-free VM equilibria. 36 So far only linear force- 
free ID VM equilibria are known— ££.£2.&LM. anc l the pseudo-particle analogy can also give 
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some insight into the types of distribution functions permitting a linear force-free solution. 

The paper is structured as follows. In Sect. [Til we present the basic general theory 
of quasi-neutral ID VM equilibria with three constants of motion, rederiving the basic 
equations given first by Mynick et al.— directly from the definitions of the basic quantities. 
In Sect. IIHI we discuss general properties of ID force-free VM equilibria based on the one- 
to-one correspondence of the mathematical problem with the motion of a pseudo-particle in 
a conservative 2D pseudo-potential and we present our conclusions in Sect. IIVI 



II. BASIC THEORY 



We assume that all quantities depend only on z and that the magnetic field has compo- 
nents B x and B y . The magnetic field components are written in terms of a vector potential 
A = (A x , A y , 0) where 

and the electric field is the gradient of an electric potential <fi such that 

E = -V<j) = -^e z . (3) 
dz 

In this case B and E automatically satisfy the homogenous steady-state Maxwell equations 
V ■ B = and V x E = 0. 

Due to the symmetries of the system (time independence and spatial independence of x 
and y) the three obvious constants of motion for each particle species are the Hamiltonian 
or particle energy for each species s, 

H s = \m s (v 2 x + vl + v 2 z ) + q s <P, (4) 

the canonical momentum in the x-direction, p xs , 

Pxs = m s v x + q s A x , (5) 
and the canonical momentum in the ^-direction, p ys , 

P ys = m s v y + q s A y , (6) 
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where m s and q s are the mass and charge of each species. 

All positive functions / s satisfying the appropriate conditions for existence of the velocity 
moments and depending only on the constants of motion, 

fa = fs(H s ,p xs ,p ys ), (7) 

solve the steady-state Vlasov equation 

v|4(E + vxB).| = 0. (8) 
or m s av 

We mention that if the same combination of values for the constants of motion allows particle 
trajectories in several distinct regions of phase space then it is in principle possible to assign 
different values to J* in each region^ 1 ^ (for an example for 2D rotationally symmetric VM 



equilibria see Ref. |42| ) . We will, however, not consider this possibility here. 

To calculate ID VM equilibria we have to solve the remaining inhomogeneous steady- 
state Maxwell equations, V • E = cr/e and V x B = yU j- Using equations (JTJ), ([2]) and (J3J) 
these equations reduce to 

d 2 4> 1 / , , ,v 
~d? = -^'A/^), (9) 

d 2 A 

— dz^ = VoJx(A x ,A y ,(f)), (10) 
d 2 A 

= W»(Ad,A/»> ( n ) 

where the source terms are the electric charge density and the current densities in the x and 
y directions, which are defined as velocity moments of the equilibrium distribution functions, 
f s , in the following way 

v-^ f°° m v 2 

a(A X} Ay,(j)) = Z^Qs fs(—k h q s <j) } m s v x + q s A X} m s v y + q s A y )d 3 v, (12) 

J-oc * 

v— \ f°° m v 2 

j x (A x ,Ay,(/)) = Z^Qs v x f s (—^- h q s (j) } m s v x + q s A X} m s Vy + q s A y )d 3 v, (13) 

J-oc * 

v-^ f°° m v 2 

j y (A x , A y , 0) = 2J g s / Vy f s (—t 1- q s( f) } m s v x + q s A x , m s v y + q s A y )d 3 v . (14) 

J-oc * 

Here, the dependence of the charge and current densities on the electric and vector potentials 
has been made visible explicitly. 
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Independently of the distribution functions f s , one can show (see Appendix [A} that the 
charge and current density components always satisfy the equations 

a<T +lf = 0, (15) 



dA 



X 



(> ° • % 0. (16) 



dAy 



il±_!MlL (,. (17) 

dAy dA x 1 ' 

These equations are completely analogous to the relation derived by Schindler and co- 
workers^ 1 ^ for the case of distribution functions depending only on the Hamiltonian and 
a single canonical momentum. 

As in that case Eqs. (j!5p - (1171) imply the existence of a potential function P, where P 
is given by 

POO 

P(A X , A y , <P) = J2 rn s v 2 J s d 3 v, (18) 



s 



the P zz component of the pressure tensor. The charge density, cr, and the current densities, 
j x and j y , are given by partial derivatives of the pressure tensor with respect to the electric 
potential, <fi and the components of the vector potential, A x and A y , in the following way 
(for details see Appendix lAl: 

OP , , 

* = -*■ (19) 

dP 

*=sc- (20) 

dP , , 

] *=8A t (21) 

Equations (fl9l) - (l2Til are equivalent to the force balance of the VM equilibrium. 24 

Throughout this work we will always be assuming a quasi-neutral plasma consisting of two 

species (electrons and ions) of opposite charge (\q s \ = e). The assumption of quasineutrality 

is justified as long as typical length scales of variation of the plasma, L, are much larger 

than the Debye length, Ad- This condition is usually satisfied very well in the majority of 

plasma systems. It should be remarked that quasineutrality does not generally imply that 

the electric field vanishes.- 14 ^ 14 ^ 

Formally, the quasineutrality condition can be regarded as the lowest order term in an 

expansion of Poisson's equation for the electric potential in the small parameter e = Xp/L. 
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In the lowest order of this expansion, quasineutrality corresponds to the condition 

dP 

<r{A x ,A y ,<f>) = —— = 0. (22) 
This equation implicitly defines the quasineutral electric potential <fi qn where 

4>qn = 4>qn(A x ,A y ). (23) 

Our potential P is then replaced by its quasineutral version P gn , where P qn is a function 
only of A x and A y \ 

Pqn(A x , Ay) = P(A X , Ay, <fiq n (A x , Ay)). (24) 

It should be noted that due to Eq. (j22]) . we have 

dPqA (dP\ + (d^A (d_P\ __(dP_\ 
dA x J Ay \dAj Ay)<p \dA x ) Ay \d^ ) AxAy \dA x ) A J 1 ) 

where it is understood that on the right hand side we evaluate P using the quasineutral 
electric potential and the subscripts indicate explicitly which quantities are being kept con- 
stant during the differentiation. An analogous equations holds for the derivative of P qn with 
respect to A y . Thus in Eqs. (1201) and (12T1) we can replace the derivatives of P by derivatives 



of P qn in the quasineutral case.— 

Using Eqs. (120]) and (T2~Tj) we can write Ampere's law as 

d 2 A x dPq n 



dz 2 ^ dA x 
d 2 A, dp 



(26) 



to^P- (27) 



dz 2 dA y 

The task of finding ID quasineutral VM equilibria has been reduced to solving these two 
coupled second order ordinary differential equations. 

It has been noticed before for particular distribution functions that the charge and current 
densities can be written as partial derivatives of a single function of the electric and the 
vector potential^^ 1 ^ as shown in Eqs. f[T9]) - (|2Tj) . However, as far as we are aware it was 
not noticed before that this is a general property of the ID VM equilibrium problem and 
that the generating function is identical to P zz . 

Grac— seems to have been the first to notice the special role played by P zz for the ID 
VM equilibrium problem, but he only investigated the special case of vanishing electric 
potential and only one non- vanishing component of the magnetic field and vector potential. 



In that case there is no need to invoke the quasi-neutrality condition and only one of the 
two equations fl20|) and fl2Tl) is non-trivial. 

Bertotti 7 generalized this to include the electric potential and one component of the vector 
potential and showed that the partial derivatives of P zz give rise to the charge and current 
density, respectively. He also discussed the assumption of quasi-neutrality for this case. 

For the = case and a set a special assumptions for the distribution functions, Lerche^ 
was able to relate the derivatives of P zz with respect to A x and A y to the respective current 
densities. For = and similar assumption for the distribution functions, other authors^ 1 ^ 
also noticed that the current density is related to derivatives of a single function with respect 
to A x and A y , but did not relate that function to P zz . 

Channel!^ showed that for = and distribution functions of the type 

f a = fsoexp(-(3 s H s )g s (p xsl p ys ), 

the current density components are given by the equations, 

1 1 \ dN 

J e + JJdA x ' 
1 1 \ dN 
J e + jJdA~ y - 



Jx = 

Jy = 



Here the plasma has been assumed to consist of two particle species (electron and ions) , (3 S = 
1/(/cbT s ) and N(A x ,A y ) is the particle density of one of the species. Another assumption 
made by Channel!^ is that N e (A x , A y ) is the same function of A x and A y as Ni(A x , A y ), 
which has implications for the distribution functions as well. It is easy to see that under 
these assumptions we obtain 

which is, of course, consistent with Eqs. (!20|) and (12T|) . 

Using force balance as an argument to show the validity of Eqs. ( Tl9l) - ( I2T1) . Mynick 
et al.— then ultimately derived the general theory presented above. Their approach is 
of course completely equivalent to ours, but we showed these equations starting from the 
velocity moments and Eqs. f|T5|) - ffT6l) . 
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Attico and Pegorarc 29 also noticed the connection between the current density and the 
partial derivatives of P zz , again for the case = and for basically the same class of distri- 
bution functions used by Channell.— In their case the distribution functions are constructed 
by linear superpositions of distribution functions of the Harris sheets type 

/, = /,„ axpHVT.) / *K„ «».) expIAt^, + u„p„)]du^u y ,. 

The authors point out that for their case Eq. ( ITTjl is a necessary condition for the existence 
of a single function from which the current density can be derived by differentiation, and 
show that it is satisfied for their class of distribution functions and that this function is P zz 
(modulo constant factors). 

One can see immediately that Eqs. (126!) and (l2"Tj) are equivalent to the equations of motion 
of a pseudo-particle with coordinates A x , A y moving in a conservative pseudo-potential 
fi Q P qn (A x , A y ). The total pseudo-energy of the pseudo-particle, E, is given by 



E 



1 fdA 
2 



dz 



+ 



1 fdA, 



dz 



+ fi P qn (A x ,A y ) = P T . 



(28) 



which, on the one hand, is also the Hamiltonian for the particle and, on the other hand, is 
the equilibrium condition for the ID VM equilibrium stating that the total pressure, P?, for 
any ID VM equilibrium is a constant. Knowledge of P qn as function of A x and A y allows us 
to predict the nature of the solution using pseudo-energy conservation without solving Eqs. 
and (j27p explicitly. Again, one or several of these properties have been noticed by a 



large number of authors for special cases (see e.g. Refs. 
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III. CONDITIONS FOR FORCE-FREE ID VM EQUILIBRIA 

The pseudo-particle analogy can be used to specify the necessary conditions that all 
force-free solutions of the ID VM equations have to satisfy. So far only linear force-free ID 
VM solutions are known^ 2 -! 2 ^^^ and the analogy could help to answer the question (see 



e.g. Refj36|) whether other, in particular non-linear, force-free solutions exist and how to 
find them. 

The force-free condition 

j x B = 0, (29) 



implies that the current density j = V x B//io is parallel to B. In the ID situation we 
discuss in this paper Eq. fl29l) can be written as 

dz \2^J ~ ^ ' 

Because the total pressure is always constant for ID VM equilibria, this implies that for a 
force-free solution P zz must also be constant. 

At first sight, this may seem to be at variance with Eqs. fl20l and (|2ip which clearly 
imply that the current density is only non-zero if the partial derivatives of P zz with respect 
to A x and A y are non-zero. A closer look, however, reveals that the condition P zz = constant 
only implies 

dP Z z _ dA x dP zz dA y dP zz _ ^ (31) 
dz dz dA x dz dA y 

for the force- free solution only. Equation ( 13TI) . which is, of course, equivalent to Eq. ( |29f) . 
can be satisfied for one solution even if the partial derivatives of P zz are non-zero. When 
translated into the pseudo-particle picture, we see that to obtain a pseudo-particle trajectory 
corresponding to a force- free magnetic field, we need a pseudo-potential (P zz ) which has at 
least one equipotential line (contour) that is also a particle trajectory. This is a necessary 
condition for the existence of a ID force- free VM equilibrium. 

Of course, finding such a potential (pressure) is not yet the complete solution of the 
problem, but only a first step. For a complete solution we also have to find the distribution 
functions giving rise to the pressure function (pseudo-potential). In simple cases this may 
be achievable by using the transform methods presented by Channel!^ and Attico and 
Pegoraro.— In more general cases, numerical methods as described by Mynick et al.— could 
be used. 

There is, however, a well-known family of pseudo-potentials that satisfies the condition 
of allowing trajectories which are identical to contours of the pseudo-potential. These are 
attractive central potentials. We have to restrict ourselves to nonsingular pseudo-potentials 
because the equivalent pressure must be positive and nonsingular. This rules out, for exam- 
ple, all potentials which are negative powers of the radial coordinate. For central pseudo- 
potentials we have not only pseudo-energy conservation, but also pseudo-angular momentum 
conservation. The pseudo-angular momentum is given by 



dAy dA x 

l-'pseudo dz dz (A'X^X AyDyj ■ 
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and is equal to the negative of the magnetic helicity density. 

An example of a distribution function resulting in a central attractive potential is given 

by : 



f s = ^exp(-(3 s H s ), (32) 



with 



V th,s 



H s = H s + ^(p 2 xs +p 2 ys ), (33) 



where a s is a dimensionless constant, which can be related to the temperature anisotropy 



of the distribution function (see e.g. Ref. |40|). We remark that the velocity space integral 
( 1181) defining P only exists if a s > —1/2. Here rio s is a constant normalizing particle density, 
P s = l/ksTg is the inverse temperature and Vth, s = (n^sPs)' 1 ^ 2 is the thermal velocity. 
For the distribution function ( 1321) the zz-component of the pressure tensor is given by 

Pzz = Ys7T exp{-(3 s q s <f))N s {A x , A y ), (34) 

a Ps 

where 

N S (A X , A y ) = n 0s exp[-r s (A 2 x + A 2 y )}. (35) 
In Eq. (1351) we have used the definitions 



n 0s = V8t^(1 + 2a s ) 1 n 0s , (36) 
r s = (3 s a s q 2 /[m s {l + 2a s )}. (37) 

The charge density is calculated using Eq. ffTUl) : 

a = ^q s exp(-(3 s q s (j))N s (A x , Ay), (38) 

s 

and the quasi-neutrality condition a = then gives 

^Tip^wAw)- < 39 > 

One can see immediately that the quasi-neutral electric field will only vanish for a choice of 
parameters such that N e (A x , A y ) cx Ni(A x , A y ). 
The quasi-neutral P zz is given by 

= P eM-r qn (Al + Al)}, (40) 
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with 



_ fie+Jk ft/(/3 e +ft) /3 e /(p e +Pi) {AU 
r — p p "'Oe "-Oi ' \^ L ) 



_ fie fit 2 
^ ' fieTfil 6 



a, 



(42) 



m e (l + 2a e ) rn.j(l + 2aj 

The x- and y-components of the current density can be calculated from Eqs. (}2T)j) and 
( !2T|) . resulting in 

j x = -2P r-^A B ex P [-r 5n (^ + ^;)], (43) 

j y = -2P r qn A y exp{-r qn {A 2 x + A 2 y )]. (44) 

At this point it is convenient to normalize all quantities. We normalize the magnetic field 
to a typical value B Q and the coordinate 2 to a typical length scale L. The components of 
the vector potential are then normalized by B$L and the normalized coefficient is defined 
by f qn = £>QL 2 r 1( j n . The pressure P is normalized by Bq/ fi and, finally, the current density 
components are normalized by HqL/B . From now on we assume that all quantities have 
been normalized in the way just described and we suppress the notation for normalization 
in what follows. 

For the case r qn < 0, P zz ,qn represents an attractive central pseudo-potential and therefore 
Ampere's law, given by the two coupled differential equations 

= 2P \r qn \A x exp[\r qn \(A 2 x + A% (45) 

= 2P \r qn \A y exp[\r qn \(Al + Al)}, (46) 
allows solutions of the form^^ 

A x = k sin az, (47) 

A y = k cos az. (48) 

These solutions represent circles in the A x -A y -p\a,ne and therefore have the property A\ + 
A 2 = k 2 = constant. This renders Eqs. ( )45l) and f j4*6i) linear because the exponential factor 
is then constant. For consistency we need that 



a = j2\r qn \P exp(\r qn \k 2 /2). (49) 



12 




FIG. 1: Left panel: A surface plot of the P zz for the case when r qn < 0. The linear force- free 
solution with A 2 + A 2 = 1 is shown above the surface plot. Right panel: The magnetic field 
components as functions of z. The magnetic field is normalised to 1 and a = 1. 

It is easy to see that 

B x = aA x = ka sin az, (50) 
By = aA y = ka cos az, (51) 

so B 2 + By = k 2 a 2 = constant as is required, and that 

j = «B. (52) 

A surface plot of P zz and the corresponding solution are shown in Fig. [IJ It is obvious that 
the linear force-free solution is a special solution of the differential equations (H5l) and (|46|) 
and that other initial conditions will produce non-force-free magnetic fields. 

Another distribution function giving rise to the same magnetic field solution has been 
presented by Channel!^ (see his case C). In this case a complete solution is available, which 
has not been discussed in the context of force-free fields, but it is obvious that for the correct 
choice of initial conditions and parameters the same linear force-free field as given above 
results. In this case P zz (pseudo-potential) has the form of a 2D harmonic oscillator potential 

P ZZ (A X , A y ) = P 00 + l -P m {A 2 x + A 2 y ), (53) 

with P o an d Poi both positive parameters. In this case Ampere's law is always linear and 
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has the form 



~p = P^, (54) 

The general solution as given by Channell is 

A x = A x0 sm(^/Pw_z + S x ), (56) 
A y = A y0 sm(>/P^z + 5 y ) . (57) 

One can see that by choosing A x q = A y o = k, 5 X = and 5 y = it/2 we recover the solution 
from above with a = V-fbi- The corresponding distribution functions are of the forrn^ 

fs(H s ,p xs ,p ys ) = exp(-(3 s H s )[f 0s + fis(p 2 xs + P 2 ys )}- (58) 

The same type of P zz and corresponding distribution function has also been found by Attico 
and Pegoraro,— but without giving explicit solutions for the magnetic field. 
It is actually straightforward to see that all distribution functions of the type 

f s = f s (H s ,p 2 s ), (59) 

with p 2 s = p 2 xs + p 2 s lead to a P zz corresponding to a central pseudo-potential. Defining 
Pxs = P cos 9, p ys = psin9, v x = vcos9, v y = -usin^, A x = A cos 9, A y = Asin9, so that 
v 2 = vl + V 2 and A 2 = A 2 X + A 2 . It is obvious that 

p 2 = m 2 v 2 + 2m s q s Av + q 2 A 2 , 

which does not depend on the angle 9. Since H s does not depend upon the vector potential, 
the integrals defining P(A X , A y ) (see Eq. (TT51) ) define a function depending only upon the 
magnitude of the vector potential, A, but not upon its direction. The corresponding pseudo- 
potential will therefore be a central potential and if it is attractive it will admit circular 
orbits, i.e. linear force-free solutions of the same type as discussed before, so obviously there 
are many distribution functions leading to the same magnetic field solution. 

Another property which is common to all force-free ID VM solutions is the following. 
Assume that a P ZZ (A X , A y ) admitting a force-free solution A x jf(z) and A y jf(z) is known, 
and that the constant value of P zz for the force free solution is Pff. Then any (positive) 
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function F(x) can be used to construct a new P ZZ (A X , A y ) admitting exactly the same force- 
free solution, by letting 

Pzz{Ac, Ay) = - ¥ ^F{P ZZ {A X , Ay)), (60) 

where F'(x) is the derivative of F with respect to its argument. Using this definition 
Ampere's law for the new P zz has the form 

= (61) 
- i4^ np ' M *' A * ))a k- (62) 

For the force- free solution we have P zz — Pff and thus Eqs. ( IBTl) and (|62l) reduce to the 
equations generated by P zz , but only for (A x , A y ) = (A x jf, A y jf). 



IV. SUMMARY AND CONCLUSIONS 



Previous work 24 on quasi-neutral ID VM equilibria has shown that generally the quasi- 
neutral ID VM equilibrium problem is equivalent to the equations of motion of a pseudo- 
particle in a two-dimensional conservative pseudo-potential and that this pseudo-potential 
is given by one of the diagonal components of the plasma pressure tensor (P zz in the present 
paper). Furthermore, the partial derivatives of P zz with respect to the electric potential 
and the two components of the vector potential generate the charge density and the two 
non- vanishing components of the current density. The equivalence of the ID VM problem to 
that of a pseudo-particle moving in a 2D potential is useful because, at least for sufficiently 
simple potentials, the general nature of the solutions can already be found from a discussion 
of the pseudo-particle problem using (pseudo-) energy conservation without explicitly solving 
the equations. 

In the present paper we have first re-derived this general theory directly using the velocity 
moments and have then applied the pseudo-particle analogy to the interesting problem of 
force-free ID VM solutions, i.e. equilibria where the current density is everywhere parallel 
or anti-parallel to the magnetic field. As we were able to show, force-free ID VM equilibria 
exist if the pseudo-potential of the equivalent pseudo-particle problem allows a particle 
trajectory which is identical to a contour of the pseudo-potential. One well-known class of 
pseudo-potentials allowing this are attractive central potentials, which have circular contours 
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and allow for circular orbits. For these circular orbits, the equivalent ID VM equilibrium 
problem is linear in all cases and leads to linear force-free equilibria. We were also able to 
show that P zz functions corresponding to central potentials can be generated by a whole 
class of distribution functions with specific properties. To the best of our knowledge, so far 
only linear force-free VM equilibria are known^^^^^^ and it is an interesting question 
whether nonlinear force-free ID VM equilibria can be found. We hope that the conditions 
formulated in this paper will help to find nonlinear force-free VM equilibria. 

One motivation for looking for nonlinear force-free ID VM equilibria is the fact that it 
is current practice to use pressure balanced VM equilibria with a guide field, mostly the 
Harris sheet£, to mimic force-free VM e quil ibria in, for example, particle-in-cell simulations 



of magnetic reconnection (see e.g. Refs. 



48 : 49). The guide-field equilibria have, 



however, the same pressure and density gradients as the neutral sheet equilibrium whereas 
proper force-free equilibria have constant density and pressure. Furthermore, a constant 
guide field does not add any free energy to the system. Also, unlike for force-free fields, 
the current density is completely independent of the strength of the added guide field. It 
would therefore be very interesting to investigate whether there are any differences between 
the reconnection process for guide field systems and for nonlinear force-free systems, for 



example how micro-instabilities found for systems wit 



49 



50 



51 



52 



53 



54 



hp 



55 



56 



asma density gradients change for 



force-free equilibria (see e.g. Refs 

The linear force-free equilibrium presented in Section II II I has been investigated using 
linear— i22iiQi£L and nonlinear— theory, but a detailed comparison of, for example, the 

magnetic reconnection process with magnetic reconnection in a corresponding neutral sheet 
has yet to be undertaken. We remark that the doubly periodic nature of this equilibrium 
may have an effect on the reconnection process compared to, for example, the Harris neutral 
sheet. It would thus be desirable to find a force- free VM equilibrium which has properties 
similar to the Harris sheet and we believe that the theory presented in this paper is a good 
point of departure for finding such equilibria. 
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APPENDIX A: DERIVATION OF RELATIONS BETWEEN THE zz- 
COMPONENT OF THE PRESSURE TENSOR AND THE CHARGE AND CUR- 
RENT DENSITY 



Following 
noticing that 



371 the relation f|T5l) can be derived directly by use of the chain rule and by 

df s 



d 6 v = 0, 



(Al) 



Vy,V z 



using integration by parts. This integral vanishes because f s has to vanish for large |v|. Also 



dfs 



m. 



dfs_ 
x V dH. 



+ 



dp x 



Pxs iPys ^ 1 ^ ' H s iPys 

using H s = m s v 2 /2 + q s (p and p xs = m s v x + q s A x . We obtain 



(A2) 



d(J , dj x / d 

+ 7u =2^ q - 



OA 



OA, 



f s d 3 v + ^- I v x f s d 6 v 



E« 



dp xs/ Hs Py 
d fs\ J3 



dh 
dH s 



Pxs iPys 



d 3 v 



8V n 



d 6 v = 0, 



which leads to Eq. ffT5]) . Equation ffT6]) can be derived in exactly the same way by replacing 
v x with v y and A x by A y . 

Equation ( TlTI) can be verified by direct differentiation of Eqs. ( 1T31) and ( fiHl) . We see that 

djx sr^ 2 f fdf s 



dA, 



s 

E 

s 

-E« 



ft 2 



1 / 



v x v v 



dv„ , 

y f v x ,v 

dfs \ ,,, 



dfs 



Pxs iPys 



d 3 v 



dH x 



d d v, 



pxs^pys 



where again the chain rule has been used in the first step and integration by parts in the 
second step. The first term vanishes during integration by parts as any admissible f s has 
to go to zero as |v| — > oo faster than any power of v. Replacing A y by A x , p ys by p xs and 
exchanging v x and v y we find with a similar calculation that 

dj y sr^ 2 f ( dfs 



dA x 



E 



Is I V x V y t 



d 3 v , 



PxsiPys 
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which shows the general validity of Eq. (TITl) . 

Equation ffl9|) can also be directly derived by differentiation of the pressure tensor in the 
following way: 



dP 



E 



m s q s / v 



d 6 v 



^ qs I Vz {^) d 



Pxs iPys 

3„ 



V 



V = — <T, 



where integration by parts has been used in the last step and 

df s 



dfs 

dv z 



m,v z 



(A3) 



Pxs tPys 



in the previous step. 

Equations (f20|) and (I2TT) are similar in structure and thus we show the derivation explicitly 
only for one of them. Differentiating the zz-component of the pressure tensor with respect 
to A x we get 



dP 



m Q / v 



dfs 



We can use Eq. ( 1A2I) to express the partial derivative of f s as 



d 3 v. 



df s 



dp x 



XS / H S ,Py S 



dv x 



m s v x 



and Eq. (1A3I) to replace m s v z (df s /dH s ) p ^ sp ^ s to obtain 



dP 

dA x 



Vy,V Z 



df^ 



dv z 



Pxs iPys 



V X ,Vy 



d 3 v 



£,9. J vj.d>v = j x . 



We have used integration by parts again in the final step and have retained only the non- 
vanishing terms. Equation (j2ip follows in exactly the same way by replacing v x by v y and 
Ac by A y . 
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